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Abstract 


We derive the decay widths for the leptonic decays of heavy charged pseudoscalars to massive 
sterile neutrinos, M± —> (A + N, within the frameworks involving the Standard Model and two- 
Higgs doublets (type II). We then apply the result to B± —> r± + “missing momentum” of the 
Belle/BABAR experimental results, in order to measure directly the relevant parameter space, 
including the mixing parameter U T n- 
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I. INTRODUCTION 


The purely leptonic decay^j] B + —> £ + u have been of great interest as a probe for new 
physics beyond the Standard Model (SM), because in the SM the decay rate can be calculated 
very precisely and new physics effects, for instance, charged Higgs contributions [I] in the 
two-Higgs doublet models [2] may appear in the tree-level contribution. In the SM, the 
decay rates are proportional to mj, the square of the corresponding charged lepton mass; 
therefore the B ± decays to e ± u and final states are highly suppressed in comparison 
to B ± —y Since the first evidence for B ± —tr^v decays was obtained by the Belle 

experiment [ L 3j, its branching fraction has been measured by Belle and BaBar [1], resulting 
in the world-average value |5] 

B exp (B ± ->• t±v) = (1.14 ± 0.27) x 1(T 4 . (1) 

This value is consistent with the SM-based prediction, 

£ SM (B ± -A t±v) = (0.7581°;°*°) X 10- 4 , (2) 

which is obtained by fitting the Cabibbo-Kobayashi-Maskawa (CKM) unitarity constraints [6j, 
at the level of approximately 1.5cr. This implies that if the measurement is improved in 
future H-factory experiments such as Belle II [7J, the comparison can clarify whether new 
physics scenarios are needed. 

Heavy sterile neutral particles (also known as “heavy neutrinos”), with suppressed mixing 
with the sub-eV SM neutrinos, appear in several new physics scenarios, such as the original 
seesaw |8] with very heavy neutrinos, seesaw with neutrinos of mass 0.1-1 TeV (HJ, or even 
scenarios with neutrino mass ~ 1 GeV |10lH2] . We will consider the reaction B ± —> t±N, 
where N stands for the heavy neutrino, any heavy sterile neutrino of Dirac type or Majorana 
type, in interpreting the measured branching fraction in the new physics perspective. Even 
if N is invisible in the detector, we can still separate B ± —y £±N signals from B + —y £ + u 
for l = e or //, because of the two-body nature of the decay whereby the momentum of the 
charged lepton in the B meson rest frame is nearly monochromatic and depends on the mass 
of N. The situation is complicated for B ± —>• t±N because there is more than one neutrino 
in the final state due to the fast decay of r ± , and the decay signature of B ± —* t±N 
becomes almost indistinguishable from the ordinary B ± —> r ± n. Therefore, we may not 
exclude the possibility that the experimentally observed signal of B ± —y may actually 

contain contributions from B ± —y t±N. 

Since the measurement of B ± —> t^u is not only an important CKM unitarity test of the 
SM, but also a very effective probe into new physics models regarding the charged Higgs, it 
is important to identify and study any unknown decay modes that can affect the measured 
branching fraction of B ± —)• as much as we can. In this paper we analyze B ± —y t±N 

both in the SM framework with a minimal extension to include N and in the frameworks 
involving two-Higgs doublets (type II). 


1 Throughout this paper, charge-conjugate modes are implied as well unless stated otherwise. 
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II. DECAY WIDTHS OF B ± -> t ± + “MISSING MOMENTUM” AND DETERMI¬ 
NATION OF RELEVANT PARAMETERS 


Massive neutrinos may be the final state particles of the leptonic decays of the heavy 
mesons (such as B ± ), if such neutrinos mix with the standard flavor neutrinos. If the 
mixing coefficient for the heavy mass eigenstate N with the standard flavor neutrino zy> 
(£ = e,n,T ) is denoted as LAvJ^] then the standard flavor neutrino ug {£ = e,fi,T ) can be 
represented as 


v e — U( Mk Uy. + U ( 

k =1 


£N 


N 


( 3 ) 


where z/*, (k — 1, 2, 3) denote the light mass eigenstates. In our simplified notation above, we 
assumed only one additional massive sterile neutrino N. The unitary extended Pontekorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix U [IB] is in this case a 4 x 4 matrix. However, our 
formulas, to be derived in this section, will be applicable also to more extended scenarios as 
well (with more than one additional massive neutrinos Nj ) . 

The decay B + —> r + N then proceeds via exchange of an (off-shell) W + SM gauge boson. 
In addition, if the Higgs structure involves two-Higgs doublets, the exchange of the charged 
Higgs H + also contributes, cf. Figs. [5] (a) and (b), as shown in Appendix A. Straightforward 
calculation, given in Appendix |Aj gives us an expression for the decay width r(F + —» r + iV); 
cf. Eq. (Al2b) in conjunction with Eqs. (A9) and (A13). 

In the decay B + —y t + u Ti within the SM with M Ut 0 and with no charged Higgs, only 


the decay mode of Fig. |5](a) contributes, and the expression (Al2b) reduces to 


Bsm (B + —» t + v t ) — —G 2 F fg\V uh \ 
o7r 


1 - 


M t 2 

Ml 


M B Ml 


( 4 ) 


Here, M B and f B are the B + meson mass and the decay constant, respectively, \V u b\ is the 
corresponding CKM matrix element, and Gp is the Fermi constant. Using this formula, 
with the values f B = 0.1906 GeV (i.e., the central value of f B = 0.1906 ± 0.0047 GeV 
Ref. [5|), M b = 5.279 GeV and t b = 1.638 x 10~ 12 s [5j, the SM branching ratio values 
Eq. (|2| from the CKM unitarity constraints J6J would then imply for \V u b\ the values \V u b\ = 
(3.44tg;}4) x 10 -3 . It is interesting that these values are very close to the values \\ 4&| = 
(3.23±0.31) x 10' 3 obtained from exclusive decays B + —y n£ + is while the inclusive charmless 
decays give significantly different values \V u b\ = (4.41 ± 0.22) x 10~ 3 [5j. 

If the above decay, with M Ut « 0, is considered within the two-Higgs doublet model type 
II [2HDM(II), | 2 ]], both modes Fig. |5](a) and (b) contribute, and the expression ( |Al2b[ ) 
reduces to 


T2HDM(II)(R + —> t + v t ) — — G 2 F f B \V u b\ 

OTT 


l - 


M 2 

Ml 


M B Mlr 2 H = r 2 H T SM (B- 


—> T V T 


( 5 ) 


2 Other notations for L7/v exist in the literature, e.g. Ur in m; Bin in ;D7iI>|- 
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where the factor r# is given in Appendix [A] in Eq. (A9a) which, in the considered case of 
B + decay, is 


r H = ~ 1 + 


Ml 

Ml 


tan 2 /3 , 


( 6 ) 


where M H is the mass of the charged Higgs H + , and tan/3 = V 1 /V 2 = v^/vu is the ratio of 
the vacuum expectation values of the two Higgs doublets (down type and up type). 

Further, if we consider the decay to a massive neutrino N, B + —> r + iV, the expressions 
(|d]) and ([5]) get extended due to M N 7 ^ 0 and due to the mixing factor U t n of Eq. (|3]) . 
However, the expressions Q and (|5]) include, in addition to the channels with the first three 
(almost) massless neutrinos Vk (k = 1,2,3), also the spurious channel with a (nonexistent) 
massless fourth mass eigenstate N' whose mixing coefficients are equal to that of the (true) 
massive N , i.e., Um- This is due to the unitarity of the 4x4 mixing matrix U appearing in 
Eq. (|3]) which implies the relation 


E ioj 2 + M = 1. 


( 7 ) 


k =1 


Therefore, any deviation from the values Q and (J5]) due to the existence of a massive 
neutrino N will be equal to r(H + — y t + N) — T(B + —> t + N)\m n =o, where the second term 
is necessary in order to avoid double counting. We will now simply denote this difference as 
T (B + —y t + N). According to the general formula (Al2b), this is then 


Tsm (B + ->• t + N ) = — G 2 F f 2 B \V ub \ 2 \U rN \ 2 M B M 2 T 


87T 


A 1/2 ( 1, 


M 2 n Ml 




Ml' M 2 b 


M 2 B M'l 


[(Ml + M 2 n ) (Ml - Ml - Ml) + 4 MlMl] 


- 1 - 


M 2 

Ml 


( 8 a) 


Ahdm(ii)(H + —> t + N) — —G 2 F fl\V u b\ 2 \U r N\ 2 M B Ml 


A 1/2 ( 1, 


M 2 n Ml 




Ml ’ Ml 


m 2 b mi 


[( M 2 T r 2 H + M 2 N l 2 H )(Ml -Ml- Ml) - Ar H l H M 2 N M 2 T ] 


- 1 - 


M 2 

Ml 


' H 


( 8 b) 


The function A appearing here is defined in Appendix |A] in Eq. (A13), and the factor Ijj in 
Eq. (A9b) which, in the considered case of B + decay, is 


, Ml 

1h = 1 + mI 


( 9 ) 


The expressions Q and ( 8 a) should then be added in the SM case, and the expressions ([5]) 

and ( 8 b) should be added in the 2HDM(I1) case, in order to obtain the full decay widths for 
± 


B ± 


—> T 


'' missing momentum. 
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We now summarize three possible cases of numerical interest for the decays B ± —> r ± + 
“missing momentum,” all in scenarios beyond the SM: 

1. If the missing momenta are only from v T of the SM, and there is charged Higgs con¬ 
tribution in addition to the SM process, the decay width is determined by Eq. (|5]) . 
Figure [l] (left panel) shows the allowed regions [shown in dark- and pale-shaded grey 
(red color online) corresponding to ±lcr and ±2cr regions, respectively] in the parame¬ 
ter space of Mh and tan f3 in 2HDM(II). To determine the allowed regions, we compare 
Eq. (Jl]) for the experimental value with Eq. (|2]) for the SM contribution in the theory 
value. 

2. If the missing momenta are due to a sterile heavy neutrino as well as the SM tau 
neutrino, while there being no charged Higgs contribution, the decay width is obtained 
by adding Eqs. Q and ( |8a| . Figure [I] (right panel) shows the allowed regions (with the 
same color assignment as described above) in the parameter space of M ^ and \U t n\ 
assuming no contributions from charged Higgs. Again, to determine these regions we 
use Eq. ([TJ) for the experimental value and Eq. ([2]) for the SM contribution in the theory 
value. We also assume, for this figure, that the sterile heavy particle N is invisible, 
and hence does not decay inside the detector. Note that the upper bound of the 
allowed region of \U t n\ goes beyond 1, which is obviously much larger than the existing 
upper bound listed in Table [1} This is mainly because the central value of the current 
world average of B(B ± —>■ r ± z/) is signihcantly larger than the SM-based calculation 
obtained from the CKM unitarity constraints, Eqs. ([TJ) and (J2|. The upper bound 
is less restrictive at low masses Mjy, because the results for the process B ± r±+ 
“missing momentum” are indistinguishable from those of SM when Mj y —> 0. 


TABLE I. Presently known upper bound estimates (c/. Ref. [13]) for \U^\ 2 = e, /L r ) f° r -Mat ~ 1, 


M n [GeV] 

\u eN \ 2 

I^VI 2 

\u tN \ 2 

» 1.0 

10~ 7 

10~ 7 

10- 2 (HU) 

« 3.0 

10" 6 

10~ 4 

io - 4 on ]) 


One thing we note is that the bound on \ U t n\ in Table [I] has been determined indirectly 
by DELPHI [T7] from the invisible decay width of Z°, i.e. e + e - —> Z° —> Nu, with 
N — u T mixing. Since v T in this reaction was not explicitly identified, the obtained 
bound is inclusive of other types of neutrinos. On the other hand, the B± —>■ r±v 
mode, where r is identified, is directly related to r — N coupling. Therefore, any 
information on \U t n\ obtained from B ± —>■ r ± n is not influenced by any other types 
of neutrinos, which makes a clear difference from the DELPHI result. In this regard, 
even though the current bound on \U t n\ from B(B ± —>■ r ± z/) is much looser than that 
of DELPHI’S, it will be of great interest if the bound can be improved or evidence for 
nonzero contribution from N is found in the future measurements of B(B ± —> t ± i'). 
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M n (GeV/c 2 ) 


FIG. 1. The allowed regions determined from the measurement of B± —> t±v in two cases: (left) the 
allowed regions in the parameter space of M H + vs. tan (3 in 2HDM(II) assuming that the missing 
momenta are only from v T of the SM; (right) The allowed regions in the parameter space of Myv 
vs. \U t n\ assuming no contributions from charged Higgs but allowing the possible contributions 
from heavy neutrino N. The dark- and pale-shaded areas (red online) correspond to ±1 a and ±2cr 
allowed regions, respectively. 


3. If the missing momenta are via a sterile heavy neutrino as well as the SM tau neutrino, 
and, there is also charged Higgs contribution, the decay width comes from adding 
Eqs. ([5]) and (8b). We will discuss this interesting case in details later. 


Up until now, a possibly important effect of suppression due to the survival probability 
was not included. Namely, if the detector has a certain length L , the produced massive 
neutrino N could decay within the detector, producing additional particles. The elimination 
of such events from the decay width T(M + —> £ + N) introduces a suppression factor Sn = 
exp[—£/( tjv7jv)]) where t ~ L//3 N is the time of flight of N through the detector (f3 N being 
the velocity), and 7 jv = (1 — ^‘n)^ 1 ^ 2 is the time dilation (Lorentz) factor. Therefore, the 
suppression factor, with which we should multiply the decay width T(H + —y r + N), is thus 


S N = exp 


L 

TN'JnPn 


~ exp 


LT n 

7 N 


( 10 ) 


where in the last relation we used /W ~ 1 and tn = 1/Tjv [= l/T(iV — » all)], in the units 
where c— 1 — h. 

Following Appendix |Bj for the decays B ± —> N and with detectors of length L ~ 1 m, 

we obtain the following: 


• If Miv ~ 1 GeV: is significantly smaller than unity only if the value of \U t n \ 2 is 

close to its present (weakly restricted) upper bound (\U t n \ 2 ~ 10~ 2 ); Sn is close to 
unity otherwise. In our present numerical analysis, we assumed SV ~ 1. 
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• If Mn ~ 3 GeV: Sn is significantly smaller than unity if at least one of the values of 
the three mixing elements \U e ' N \ 2 is close to its present upper bound (| U ei \r | 2 ~ 1CT 6 , 
or II/^tvI 2 10 4 , or | U t n | 2 ~ 10 4 ); Sn is close to unity otherwise. In our present 
numerical analysis, we assumed Sn ~ 1. 


III. NUMERICAL ANALYSIS AND DISCUSSIONS 

In this Section, we discuss implications from possible future measurements. Although 
Fig.0 (right) shows that we can set some constraints on the heavy neutrino mass and its 
coupling through the measurement of B{B ± —> r^u), the existing uncertainty is too large. 
Since the Belle II experiment is aiming to increase the data sample by more than a factor of 
50, we may expect that the uncertainty of B(B ± —>■ r ± z/) can be reduced at least by an order 
of magnitude. Therefore, in the discussions below, we will assume that the experimental 
uncertainty is improved by a factor of 10. 

First we assume no contribution from charged Higgs such as in 2HDM(II). Consequently, 
we consider only the diagram shown in Fig. [5](a) in Appendix A where both ordinary u T and 
heavy neutral particle N contribute. Depending on the values of Uin, Mn and the detector 
size, the produced sterile neutrino N can decay within or beyond the actual detector. When 
N is produced in the decay B + —> t + N, and if it decays within the detector, the main 
signature of N will be N —)■ 1 + ti~ (if N is Majorana) or N —>■ 1~tt + (if N is Dirac or 
Majorana), which will show experimentally as a resonance in M(/ ± 7r =F ). In order to maximize 
experimental sensitivity, we should analyze both the invisible mode and the visible decay 
modes for N, and combine the corresponding signal yields. Consequently, the value of Uin 
shall be determined by the combined yields with appropriate corrections for efficiency and 
subdecay branching fractions. While the invisible mode of N will be analyzed by following 
the existing B + —> t + v t analyses of Belle and BABAR, some of the visible decay modes of N 
should be explicitly analyzed in the future measurement. Based on our expectation of the 
branching fractions of the visible modes (see Appendix C) and the survival suppression factor 
(see Appendix B), the correction factor for the undetected signal events can be obtained for 
each assumed values of Uin and m^. 

Figure [2] shows the allowed regions in the parameter space of Mn and \U t n\ in this case. 
The shaded area (red online) corresponds to ±lcr (dark) and ±2cr (pale) allowed regions. 
For the plot on the left panel, the central value of B(B ± —>■ t±v) is taken from Ref. |5j, 
i.e., Eq. ([Tj) , but with the uncertainty reduced by a factor of 10. The plot on the right 
panel assumes that the central value is equal to the value predicted by the CKM unitarity 
constraint [6], i.e., Eq. (J2|, again with tenfold reduction of uncertainty (i.e., ±0.027 x 10 -4 ). 
In both cases, the allowed regions are determined by comparing the expected experimental 
outcome to the theory value, where the SM contribution is taken from Eq. ([ 2 ]). 

From Fig. [2] (left), it is evident that we will need an additional contribution from, e.g., 
B + —> t + N if the central value of the current measurement of B(B ± —>• r ± zz) stays the same 
while a substantial reduction of the measurement uncertainty is achieved. Furthermore, 
comparison of Fig. [2] (left) with the inclusive upper bound in Table [I] implies that such 
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a scenario would require additional new physics, e.g. 2HDM(II) charged Higgs exchange 
contributions. Note that, in our numerical analysis that follows, we assume the suppression 
factor Sjy = 1. However, with \U t n\ 2 ~ 10~ 2 , 7 # = 2 and L = 0.1 m, we get Sn ~ 1 
(for M/v = 1 GeV), but Sn ~ 0 (for Mn = 3 GeV), which means that the heavy sterile 
neutrino N is likely to decay within the detector if Mn > 3 GeV. Therefore, it is important 
to consider visible decay modes of N as well as the invisible mode. 



FIG. 2. The allowed regions in the parameter space of Mn and \U t n\ assuming no contributions 
from charged Higgs, where the Tier and ±2cr allowed regions are displayed by dark and pale shades 
(red online), respectively: (left) the central value of B(B^ —y t±v) is taken as the current world 
average Eq. ([!]), while tenfold reduction of uncertainty is assumed; (right) the central value is taken 
to be the value predicted by the CKM unitarity constraint, also assuming tenfold reduction of 
experimental uncertainty. In all cases, the comparison is made to the value determined from the 
CKM unitarity fitting Eq. (J2|. 

On the other hand, if we consider the case where there is no contribution from unknown 
heavy neutrino N, we note that the parameters M H + and tan /3 can be much further con¬ 
strained if the B(B ± —>■ r ± z/) uncertainty is improved, e.g. by a factor of 10. Figure [3] shows 
the allowed regions in the parameter space of M^+ vs tan /3 of 2HDM (type II) while assum¬ 
ing no contributions from heavy neutral particle N. The ±lcr and ± 2 cr allowed regions are 
displayed by dark and pale shades (red online), respectively. The left panel plot uses, for 
the central value of B(B ± —> r^u), the current world average and assumes tenfold reduction 
of uncertainty. For the right panel plot, we consider the case of the central value being 
identical with the present value predicted by the CKM unitarity constraint and the experi¬ 
mental uncertainty is reduced by a factor of 10 compared to the current value ±0.27 x 10 -4 . 
In both cases, the comparison is made to the value determined from the CKM unitarity 
constraints Eq. ([2]). Figure [3] (left) shows that 2HDM(II), for each tan/3, has a very narrow 
interval of the corresponding allowed values of Mu = M#(/3) if the central experimental 









value of B(B ± —> t±v) remains approximately unchanged and the experimental uncertainty 
is reduced tenfold. 



FIG. 3. The allowed regions in the parameter space of Mjj+ vs. tan/3 in 2HDM (type II), assuming 
no contributions from heavy neutral particle N, where the Tier and ±2u allowed regions are dis¬ 
played by dark and pale shades (red online), respectively: (left) the central value of —» t^u) 

is taken from the current world average, while tenfold reduction of uncertainty is assumed; (right) 
the central value is taken to be the value predicted by the CKM unitarity constraint, also assuming 
tenfold reduction of experimental uncertainty. In both cases, the comparison is made to the value 
determined from the CKM unitarity constraints Eq. ([2]). 


Now, let’s consider the case where both charged Higgs H ± and heavy neutrino N con¬ 
tribute to the measurement of B ± —> t ± v 1 based on the decay rates of Eqs. ([5]) and ( |8b| ). 
Figure [4] shows a few exemplary cases. As in the cases of Fig. [2j we assume, in Fig. [4j 
that both invisible and visible decays of N are analyzed with appropriate corrections being 
applied to the signal yields to obtain the necessary branching fraction. For each of the two 
plots in the top panel, we choose a point in the parameter space of Mn vs \U t n\ and show 
the allowed region in the parameter space of M^+ and tan /3. For the bottom panel, we 
choose points in the space of M H + vs tan/3 and show the allowed region in \U t n\ vs M^. 
The two plots in the left panel correspond to the case where we choose points within the 
allowed region, while we choose points outside the allowed region for the two plots in the 
right panel. For Fig. |d](a), we choose M N = 1.0 GeV/c 2 and \U t n\ = 0.6. On the other hand, 
for Fig. |d](b) we choose Mjv = 1.0 GeV/c 2 and \U t n\ = 0.5. Although it may seem a small 
difference between the two cases, the resulting allowed regions shown in the M#+-vs-tan/3 
space is clearly different. Similarly, we choose M H + = 200 GeV/c 2 for both plots in the 
bottom panel of Fig. [4j but tan (3 = 56.5 (allowed) for the left and tan /3 = 55 (excluded) for 
the right. Again, we see clear difference between the two cases. 
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M n (GeV/c 2 ) 



FIG. 4. The allowed regions under the assumption that both H + and N contribute to the measured 
value of —» t±v, based on the decay rates of Eqs. ^ and (8b). Top left: the allowed region in 
the parameter space of M H + and tan/3 when M n = 1.0 GeV/c 2 and \U t n\ = 0.6. Top right: the 
allowed region when Mn = 1.0 GeV/c 2 and \U t n\ = 0.5. Bottom left: the allowed region in the 
parameter space of \U t n\ and Mn when M#+ = 200 GeV/c 2 and tan/3 = 56.5. Bottom right: the 
allowed region when M H + = 200 GeV/c 2 and tan/3 = 55. 


IV. SUMMARY AND PROSPECTS 

We have seen, in Figs. [l]|4j that the existing measurements of B ± —> t ± v can be rein¬ 
terpreted by including the possible contributions from a heavy neutral sterile particle N. 
However, with the current experimental uncertainty, the case is not clear yet. On the other 
hand, with the upcoming next-generation measurements from flavor physics facilities such 
as Belle II, the experimental uncertainty can be greatly reduced, as was indicated in Figs. 
2-4. In that case we can set interesting constraints on parameters of 2HDM(II) and on heavy 
neutrino N in the range Mn /$ 3 GeV/c 2 . For instance, if the current experimental average 
value of B(B ± —> ^v) stays approximately the same but the experimental uncertainties are 
greatly reduced, it can be a strong indication that we may need a charged Higgs such as in 
2HDM(II), and/or possibly a heavy neutrino N. 
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In the decays B ± —>• £ ± +N when £ = e or /i, the problems of double counting encountered 
in Sec. 0do not appear, because the kinematics allows to distinguish such decays from those 
of B ± —> + Uk (k = 1, 2, 3) as mentioned in the Introduction. The relevant formulas for 
such cases are thus 

r SM (B + -4 tN) = ^G 2 F f 2 B \V ub \ 2 \Um\ 2 \' /2 (l, 

x T- [(A/f + - Ml - Ml) + 4 MlMl ], 

ivi B 

r 2H DM(ii )(B+ £ + N) = ±G 2 F f 2 B \V ub \ 2 \U m \ 2 \ 1/2 (l, ||, ^ (11b) 

x + M 2 N l 2 H )(M 2 B -Ml-M 2 ) - 4 r H l H M 2 N M 2 } , 


Ml 

Ml" 


Ml 

Ml 


(11a) 


where £ = e or /r. Although the mass of £ is in this case small or almost zero implying the 
helicity suppression in the three-generation case, the presence of a massive neutrino M^ ~ 1 
GeV may make the decays GD appreciable, depending certainly on the mixing strength 
\U fN \ 2 . In the e + e _ i?-factory experiments where B mesons are produced via e + e~ —)■ 
T(4S') —> BB process, the presence of heavy neutrino N in the decays B + —> £ + N {£ = e, /i) 
will be distinguishable from B + —)■ £ + u^ by the momentum of £ + in the rest frame of B + . 

This study can be also extended to other decay modes such as B —> ZA*V + z/ and B + —> 
£ + u (£ = e or /i). By combining B ± —>■ r ± z/ and B —> D^t + v together, the sensitivity of 
searching for N can be even more enhanced. Moreover, it is expected that B ± —s- /Mv decays 
can be observed in the Belle II experiment. This decay, unlike B ± —>■ r±v, is a two-body 
decay mode of 5 , hence the hnal-state charged lepton (/i ± ) has a nearly monoenergetic 
distribution in the rest frame of the B meson. If a heavy neutrino N, in addition to of SM, 
also contributes to this decay, it will change the energy distribution of ^ and its effect can 
be measured experimentally. In the case of B ± —>■ e ± u, the SM expectation is very low, well 
beyond the sensitivity of Belle II. Nevertheless, if a heavy neutrino exists and contributes to 
B ± ->■ e ± z/, it can enhance the branching fraction to be within the experimental sensitivity 
of Belle II. 
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Appendix A: Derivation of the Decay Widths T(M ± —> &N) 



FIG. 5. The decay M + — > N £ + via the exchange of (a) W + and (b) H + . 


In this Appendix, we derive the decay width of the process M + —» £ + N of Figs. [5^a) 
and[5](b). The first decay mode Fig. |5](a), involves an exchange of (off-shell) W + . Direct 
calculation gives for the contribution ■ of the W + exchange to the reduced scattering 
(decay) matrix (assuming Mw Mm) 

T tw> = u; N v^M( 0 \D 7 "(i - ls )U\M+) [% 7 ,(1 - 75 )« ( ] , (Al) 

where D and U denote the two valence quarks of the pseudoscalar meson M + , Vup is 
the corresponding CKM mixing matrix element, and Gp = (1 /( 4 a / 2 )) (g / Mw) 2 = 1.166 x 
10" 5 GeV “ 2 . 

In the model with 2HDM(II), the couplings of the charged Higgses H ± with fermions are 
similar to those of [2]. The relevant parts of the Lagrange density are 

^H+qq = - ^ H+V UMk [ cot @UjMjj. (1 — 7 5 )D k + tan / 3UjM Dk (1 + 75 )D k ] 

V W j,k 

+H.c. , (A2a) 

C H ±m = |— H + U* m [cot PNM n { 1 - 75 Y + tan $NM t ( 1 + l5 )i] 

2v 2M\y £ 

+H.c. . (A2b) 


The first density is for the coupling with quarks, and the second for the the coupling with 
the sterile massive neutrino N and the three charged leptons ( l = e,/q t). Here, tan/3 = 
Vi/v 2 = Vd/vjj is the ratio of the vacuum expectation values of the two-Higgs doublets (down 
type and up type). 

In analogy with the case of W + exchange, we obtain the contribution of the H + 
exchange to the reduced scattering (decay) matrix T (assuming M H M m ) 

T {h) = u mV^ D ^(- 1)^- [(0|D(1 - r y 5 )U\M + )m D tan/3 + (0|D(1 + 75 )t/|M+)m t7 cot/3] 
x {[tfjv(l 75 )^] M n cot /3 + [iZjv(l T 75 )^] Mi tan /3} , (A3) 


where Mp is the mass of H + . 


The expressions (Al) and (A3) can be further simplified by the axial-vector current and 
pseudoscalar relations 


(0\D^(l- l5 )U\M + )=ifMP n M , 

— M 2 

(0|D(1 =f 7 5 )D|M + ) = =r/m — , 

m D 


(A4a) 

(A4b) 
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where Jm, Mm and pm are the decay constant, mass and 4-momentum of M + , respectively. 
In Eq. (A4b), the mass of mu was neglected in comparison with mu- For example, for B + 
we have D = b and U = u. 


Using the relation (A4a) in the expression (Al), and using the relations 

j4r = j/jv + l/i , 

u N rf N = M n u n , pfpvt = -Meve , 


(A5a) 

(A5b) 


we obtain the following form for the reduced scattering (decay) matrix element for the 
IU + -mediated decay: 


T (W> = if M U; N VJ D ^| {+M n Ml - 75H - Ml [B„(l+ 7sM} 


(A6) 


Using the relation (A4b) in the expression (A3), and neglecting there the term propor¬ 
tional to mu cot/3 (<§; mo tan/3), the reduced scattering (decay) matrix element for the 
i/ + -mediated decay becomes 


T m = i/»C/; N VJ D %^Atan/3 


s /2 Ml 

x {[u N (l - 75 )^] M n cot f3 + [mjv(1 + 75 ) 77 ] M e tan/3} . 


(AT) 


Combining Eqs. (JA6J) and (A7), we obtain finally the reduced scattering (decay) matrix 
element T for the decay M + —y N£ + 

7 - _ q-{W) _j_ ‘jiH) 

(A8a) 


— ifMU% N VuD^J= {[^iv(l + 75 )] M ( r H + [u N (l — 75 )] M N l H } , 


where the coefficients are 


r H = -1 + tan 2 /3 , 

/ + 
1h ~ +1+ mI■ 


(A9a) 

(A9b) 


The square of this matrix element (summed over the helicities of the two final particles) 
then gives 


(|T| 2 ) = ±G 2 F f 2 M \U m \ 2 \V UD \ 2 [(. M 2 r 2 H + M 2 N l 2 H )(p N ■ p e ) - 2 r H l H M 2 M 2 N \ , 
where we have 


(PN-Pl) = g ( M J 


M M 2 r M 2 ) 


In the rest frame of M + we then have for the decay width 


T(M+ ->• £ + N) = 


-A 1 / 2 ( 1, 


M 2 n M 2 


16t tM m V M 2 m 


(\T\ 2 ) 


(A10) 

(AH) 

(Al2a) 


— -^G F f M \Um\ \Vud\ 


2 A x / 2 ( 1, 


M 2 M 2 


I’M 2 ’ M 2 

V 1V1 M 1V1 M 


x— [(M 2 r| + M 2 n 1 2 h ){M 2 m - M 2 n - Ml) - 4 r H l H M 2 N M 2 ] , (Al2b) 

1V1; ' 


M 
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where we used the notation 


Hvi.m, vs) = X + vl + vl - ' h htn - - 'hhih ■ 


(A13) 


The expression (Al2b), in conjunction with the expressions (A9) and (A13), is the explicit 
expression for the decay width T(M + — > N£ + ) in the rest frame of M + , in terms of the 
masses Mm, Mn, Mi, Mh and tan/3 = vjj/vd- 

It is straightforward to check that in the case of Mn = 0 and \Uin\ — 1, the obtained 
formula (Al2b) reduces to the formula obtained in Ref. |T|. Further, if Mh —> oo (i.e., no 
charge Higgs interchange), the formula (A12b) reduces to Eq. (2.5) of Ref. [15]. 


Appendix B: Survival Suppression Factor 


We first mention that the factor of “nonsurvival” probability Pn = 1 — S/v has been 
discussed and investigated in the literature for the processes where the intermediate on-shell 
particle (such as N ) is assumed to decay within the detector, cf. Refs. [15] 1T81 - I24] . Here in 
this Section we follow the notations and results of Ref. [22]. The total decay width Tn of 


sterile massive neutrino N appearing in Eq. (10) can be expressed as 


where 


r N = y(m n )k , 

GlMl 


r (m n ) = 


_ AT F 1v± N 


967 t 3 


(Bl) 


(B2) 


and the factor /C is proportional to the heavy-light mixing factors \U e ' N \ 2 [where U e * N appear 
in the relation 


IC(M n ) = /C — We AT \U eN \ 2 + jVftN \U[j,n\ +W-JV \U tN | 2 . (B3) 

Here, the coefficients J\f e > N (Mat) — {( = e,/r, r) are the effective mixing coefficients. 

We have Af £ ' N ~ 10°-10 1 . They are functions of the mass M^ and were numerically evaluated 
in Ref. [22, for the Majorana neutrino N, on the basis of formulas of Ref. [25] . The numerical 
results for the Dirac N were included in Ref. [ 23] , 

In the ranges of M]y typical for the B ± —>■ N£± decays, i.e., for M^ & 1-4 GeV, we have 
Af e N ~ Af lt N ~ 8 and Af T N ~ 3, and therefore 

K. ~ 8(|T/eiv1 2 + \U/j,n\ 2 ) + 3|f/ r jv| 2 {B,B C decays) . (B4) 


The estimate (B4) is valid for Majorana N. For Dirac N it is somewhat lower, but the 
difference can be ignored at the level of precision of the estimate, for 1 GeV < Mn < 3 GeV. 

We refer to Ref. [22] (and references therein) for more details on these results. Further¬ 
more, approximate values of the presently known upper bounds for the squares of the mixing 
elements, in this range of masses Mn, are given in Table [I] (cf. Ref. [13]). 
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The survival factor (10) can be rewritten as 


Sn = exp ( — j— J = exp 




(B5) 


where Ln is the decay length, and Ln is the canonical decay length (canonical in the sense 
that it is independent of the mixing parameters U ( > N ) 


J N 


T 1 — 

J^N ~ 


— L N K, , 

r (m n ) 


In 


(B6a) 

(B6b) 




where T(Mn) is given in Eq. (B2). The inverse canonical decay length L N , for 7 n = 2, is 


-^-1 


given in Fig. p as a function of Mn- Specifically, we obtain L N « 10 2 m , 3 x 10 4 m , for 
Mn = 1 GeV, 3 GeV, respectively. Combining this result with the results (B4) and Table [TJ 
we obtain for the effective inverse decay length [appearing in the survival factor Sn of 
Eq. (B5)| the following estimates, in units of m . 

lGeV) 


T -1 
L N 


(Mi 


N 


~ 0.8 x 10 d |G eiV / + 0.8 x lO^U^l 2 (+2 x W> 2 \U tNj V) 
< 10” 4 + 10“ 4 (+ 10 °) , 


L~n (M n « 3GeV) » 3 x 10 5 1 t/ e7Vj . | 2 + 3 x lO 5 ^, | 2 (+1 x 10 5 |G TiV ,| 2 ) 


< 10 ° + 10 ° (+ 10 °) 


(B7a) 

(B7b) 



M w ( GeV) 

FIG. 6. The inverse canonical decay length L N = T(Mn)/'In, in units of m^ 1 , as a function of 
the neutrino mass Mn, for the choice 7 tv [= (1 — P‘n)~ 1 ^ 1 \ = 2. 


Appendix C: Branching ratios for semihadronic decays of neutrino 

Here we summarize some of the formulas for the decay widths and branching ratios for 
the decays of a heavy neutrino N into hadrons [T3| I25|. Comparatively appreciable channels 
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are decays into light mesons (with mass Mp < 1 GeV), which can be presudoscalar (P) or 
vector (V) mesons: 


2T(N rp + ) = \B m \ 2 ^M 3 N f 2 \V P \ 2 F P (x £ ,x P ) 

o7r 

r (N u ( P°) = \B £N \ 2 ^fM 3 N f P (l - 4) 2 

o47r 

2T(N rv + ) = \B m \ 2 ^M 3 N f 2 \Vy\ 2 F v (x £ ,xy ) 

07T 

p (N -4 v e V°) = \B m \ 2 ^M 3 N f 2 K 2 y(l - 4) 2 ( 1 + 24). 


(Cla) 

(Clb) 

(Clc) 

(Cld) 


Here, £ = e, /i, r stands generically for a charged lepton. The charged meson channels above 
were multiplied by a factor 2, because if N is Majorana neutrino both decays N —> i~M + 
and iV —> £ + M~ contribute (M = P, V ) equally. The factors fp and fy are the decay 
constants, and Vp and Vy are the CKM matrix elements involving the valence quarks of the 
corresponding mesons. 

In Eqs. (Cl), the following notation is used: Xy = My/Mpr (Y = £,P,V), and the 
expressions Fp and Fy are 


where 


F P (x,y) = A 1/2 (l,4 y 2 ) [(1 + x 2 )(l + x 2 - y 2 ) - 4x 2 ] 
F v (x, y) = A 1/2 (l, x 2 , y 2 ) [(1 - x 2 ) 2 + (1 + x 2 )y 2 - 2 y 4 \ 

A(?/i, 2/2,2/3) = 2/1 + 2/2 + 2/1 - 22/12/2 - 22/22/3 - 22/32/1 • 


(C2a) 

(C2b) 

(C3) 


The (light) mesons for which formulas (Cl) can be applied are: P ± = 7r , iG* 1 ; P° = 
C ± = p ± ,A'* ± ; H° = If M N > 1 GeV, the neutrino N can 

decay into heavier mesons, and the decay widths for such decays can be calculated by using 
duality, as decay widths into quark pairs; nonetheless, such decay modes are in general 
suppressed by kinematics and are not given here. The corresponding branching ratios are 
obtained by dividing the above decay widths by the total decay width T^ of N. 
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